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Abstract 

Using the Genefishing method, we identified seven potential regulatory genes involved in the process of scale 
morphogenesis in fishes. We further characterized a novel solute carrier protein gene {CcSLQ, from the common 
carp which is differentially expressed in mirror carp and Jianli. The ORF encodes a peptide of 298 amino acids with 
a molecular mass of 31.5 kDa and a theoretical isoelectric point of 7.49. ScanProsite analysis indicated that it is a 
putative solute carrier protein that contains a substrate binding site. CcSLC was detected in carp embryos by in situ 
hybridization in the 70%-epiboly, 6-somite, and 14-somite embryonic stages. Gene expression stopped at the long 
pec stage. However, CcSLC25a5 was re-expressed during the initiation of scale formation in the regions that were 
scale covered. These findings provide novel insights into the features of early carp embryo and scale development. 

Keywords: Cyprinus carpio Linnaeus; Scale development; Solute carrier protein; Scale initiation; Embryo 
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Introduction 

Membrane transporters are the gatekeepers for all cells 
and organelles, controlling uptake and efflux of crucial 
compounds such as sugars, nucleotides, inorganic ions, 
and drugs (Hediger et al. 2004). They are responsible for 
substrate movement across both cytoplasmic mem- 
branes of cells and internal membranes of organelles 
(Sreedharan et al. 2011). Transporters can be divided 
into ABC transporters, pumps, ion channels, water 
channels, and solute carriers. Membrane bound proteins 
represent about 27% of the entire human proteome. 
Among the membrane bound proteins, the SLC trans- 
porters are the second largest group after G protein 
coupled receptors (Lagerstrom 2008, Almen et al. 2009). 
Transporters can be divided into two families, passive 
and active transporters. The active transporters use 
diverse energy-coupling mechanisms to allow the move- 
ment of molecules across a membrane against a concen- 
tration gradient. The passive transporters, also known as 
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facilitated transporters, allow passage of solutes (e.g., 
glucose, amino acids, urea) across membranes down 
their electrochemical gradients (Hediger et al. 2004). 

Appreciation of the role that transport proteins play in 
the absorption, distribution, and elimination of a wide 
variety of drugs in clinical use is increasing. As the lar- 
gest group of secondary transporters, SLC transporters 
are becoming the focus of an increasing number of stud- 
ies because they control transmembrane movement of 
many types of important substrates. The human genome 
contains approximately 360 unique SLC protein genes 
grouped into 48 families (Ren et al. 2007; Fredriksson 
et al. 2008). Approximately 19 of the SLC gene families 
have been reported to transport xenobiotics including: 
organic anion polypeptides {SLCO), oligopeptides {SLC1S) 
(Russel et al. 2002; Brandsch et al. 2008; Dobson and Kell 
2008; Rubio and Daniel 2008), organic anion/cations 
(SLC22) (Koepsell et al. 2007; Ciariboli 2008), and organic 
cations (SLC47) (Tanihara et al. 2007; Moriyama et al. 
2008; Matsushima et al. 2009). 

The SLC2S gene encodes mitochondrial carriers (MCs), 
which are membrane-integrated proteins that localize to 
the inner membranes of mitochondria and catalyze the 
translocation of solutes across the membranes (Plamieri, 
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2004). The MCs provide a critical link between the mito- 
chondria and the cytosol by facilitating the flux of solutes 
through the permeable barrier of the inner mitochondrial 
membrane. The substrates transported by the MCs range 
from the smallest H+ to the largest ATP molecule, imply- 
ing that they have a broad array of functions in diverse 
metabolic processes. Defects in MC genes lead to several 
diseases such as type II citrullinaemia (SLC25A13; OMIM 
215700), hyperornithine-hyperammone- homocitrulline- 
mia (HHH) syndrome (SLC25A15; OMIM 238970), Stanley 
syndrome (SLC25A20; OMIM 212138), Amish microceph- 
aly (SLC25A20; OMIM 607196), and autosomal dominant 
progressive external ophthalmoplegia (adPEO) (SLC25A4; 
OMIM 157640). The complete amino acid sequence of 
the ATP/ADP carrier was identified in beef heart mito- 
chondria (Aquila et al. 1982; Aquila et al. 1985). 

Post-genomic era studies have enabled us to identify 
many more mitochondria carrier families (MCFs) simul- 
taneously without laborious cloning or purification proce- 
dures. Although much is known about the characteristics 
and functions of MCFs in human and plants, their bio- 
logical roles in fish remain unknown. In our studies, we 
cloned the CcSLC2SaS {Cyprinus Carpio SLC25a5) gene 
using Genefishing kits from the skins of the mirror carp, 
which has interspersed scales, and the Jianli, that has full 
scales. The expression pattern of SLC25a5 during different 
developmental stages was determined by whole-mount in 
situ hybridization. 

Materials and methods 

Animals 

Mirror carp and Jianli {Cyprinus carpio Linnaeus) were 
cultivated at Experimental Station of the Wuxi Fresh- 
water Center, Jiangsu, China. The mirror carp was de- 
rived by domesticating the common carp and selecting 
for a scale-reduced mutation fgfrla (Rohner et al. 2009). 
The skin tissues from mirror carp and Jianli were har- 
vested with forceps and immediately homogenized in 1 ml 
Trizol (Invitrogen). 

First-strand cDNA synthesis 

Total RNA extracted from the skin tissues using Trizol re- 
agent (Invitrogen) was used to synthesis the first-strand 
cDNA. Subsequent reverse transcription was performed 
according to the manufacture's protocol (Seegene, Seoul, 
South Korea). The final reaction volume was 20 uL and 
contained: 3 ug of purified total RNA, 4 uL of 5x reaction 
buffer, 5 uL of dNTPs (2 mM each), 2 uL of 10 uM dT- 
ACP1 (5 ' -GTCTACCAGGCATTCGCTTCATXXXXXGC 
CATCGACC-3'), 0.5 uL RNase inhibitor (40 U/uL; 
Invitrogen, USA), and 1 uL of reverse transcriptase (200 U/uL, 
Invitrogen). First-strand cDNAs were diluted using 80 uL 
of DNase-free water for GenefishingTM PCR, and stored 
at -20°C. 



ACP (Annealing Control Primer)-based Genefishing PCR 

DEGs (Differential Expressed Genes) were screened by 
ACP-based PCR methodology using the Genefishing 
DEG Kits (Seegene). Briefly, second-strand cDNA was 
synthesized at 50°C during in the first-stage PCR reac- 
tion. The final reaction was conducted in a 20 uL vol- 
ume containing: 3-5 uL of diluted first-strand cDNA, 1 
uL of dT-ACP2 (10 uM), 1 uL of 10 uM arbitrary ACP 
(Hwang et al. 2005), and 10 uL of 2x Master Mix 
(Seegene). The PCR protocol for second-strand synthesis 
was: one cycle at 94°C for 5 minutes, followed by 50°C 
for 3 minutes, and 72°C for 1 minute. Once the second- 
strand DNA synthesis was completed, a second-stage 
PCR amplification protocol was conducted that consisted 
of: 40 cycles of 94°C for 40 seconds, 65°C for 40 seconds, 
and a 5 minute final extension at 72°C. The amplified 
PCR products were separated in a 2% agarose gel and 
stained with ethidium bromide. 

Cloning and sequencing 

PCR bands indicating genes with differential expression 
were extracted from the gel using a DNA extraction kit 
(Zomanbio, China). The bands were directly cloned into 
a pEASY-T vector (Trans, China) according to the man- 
ufacturer's instructions. The cloned plasmids were 
sequenced. 

Whole-mount in situ hybridization 

RNA probes were prepared from a 206 bp CDS (Coding 
Sequence) region of the gene SLC2SaS in common carp 
and labeled with digoxigenin-UTP using T3 or T7 RNA 
polymerase (T3 for production of the antisense probe, 
T7 for the sense probe). The embryonic and develop- 
mental stages of the embryos used for whole-mount in 
situ hybridization were assessed using haf (hours after 
fertilization) and various morphological criteria (Kane 
and Kimmel 1993) as described by Westerfield (1993). 
The RNA probes were hybridized to the tissue overnight 
at 65°C. The embryos and juvenile fish from each devel- 
opmental stage were imaged using an Olympus BH-2 
microscope (Olympus Optical, Tokyo, Japan). The pri- 
mers used to create the probes were: Forward: 5 -TGGG 
TAACTGCTTGGTGAAGATCTCC-3 ' , and Reverse: 5 -A 
CCAGCAACAGCAGTCACAGTCTGA-3 ' . 

Results 

The mirror carp and Jianli have differential gene 
expression in skin tissues 

To identify the differentially expressed genes that are as- 
sociated with development of the skin and appendages, 
skin samples from mirror carp and Jianli were assayed 
by Genefishing. The Genefishing assay used an anchored 
primer in combination with 20 arbitrary primers (Arbitrary 
ACP, Annealing Control Primer) (Table 1). We obtained 
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seven DEGs from all ACP primers (Table 2), among these 
DEGs, two PCR products from ACP28 and ACP29 primers 
were identified that had significandy different expression 
levels (Figure 1). The differentially expressed bands were 
subcloned into pEasy-T3 vector and sequenced. The se- 
quences were compared using the blast program in the 
NCBI sequence database (http://wei.sohu.com/20121011/ 
n354628269.shtml). The blast results showed that the dif- 
ferentially expressed genes identified using the primers 
ACP28 and ACP29 were homologues of the zebrafish 
SLC25a5 (So/ute carrier) and TPT1 (Tumor protein, 
translationaly-controlledl) genes, respectively. We have 
named them CcSLC25aS (The Common carp SLC25a5) 
and CcTPTl (The Common carp TPT1). 



Table 2 Seven candidate GOIS (gene of interest) by 
genefishing 



GOI number Protein 

1 Solute carrier family 25 alpha, member 5 (slc25a5) 

2 Hairless protein 

3 Cyprinus carpio translationally-controlled 

tumor protein mRNA, complete cds 

4 Danio rerio myosin heavy chain, fast skeletal 

muscle-like, 85 bp (89%) 

5 Actin, alpha, cardiac muscle 1 b [Danio rerio], 

367/389 (94%) 

6 Dictyostelium discoideum AX4 hypothetical protein, 

84/87 (97%) 

7 Fibroblast growth factor 4 



CcSLC25a5 gene structure and its isoforms 

One of the identified genes, CcSLC2Sa5, was selected for 
further analysis. We designed the primers located in the 
conserved region of the gene using the partial reference 
sequence of CcSLC25a5 isoform 1. Nine different isoforms 
of the CcSLC25a5 gene were identified in the common 
carp and Jianli. The genes encode six different SLC25a5 
protein isoforms due to codon degeneracy. Four isoforms 
were found in the common carp and two isoforms were 

Table 1 Primers used in genefishing for amplifying the 
differential expressed genes in skin of jianli and mirror 



carp 


Primers 


Primer sequence 


ACP1 


GTCTACCAGGCATTCGCTTCATXXXXXGCCATCGACC 


ACP2 


GTCTACCAGGCATTCGCTTCATXXXXXAGGCGATGCC 


ACP3 


GTCTACCAGGCATTCGCTTCATXXXXXCCGGAGGATG 


ACP4 


GTCTACCAGGCATTCGCTTCATXXXXXGCTGCTCGCG 


ACP5 


GTCTACCAGGCATTCGCTTCATXXXXXAGTGCGCTCG 


ACP6 


GTCTACCAGGCATTCGCTTCATXXXXXGGCCACATCG 


ACP7 


GTCTACCAGGCATTCGCTTCATXXXXXCTGCGGATCG 


ACP8 


GTCTACCAGGCATTCGCTTCATXXXXXGGTCACGGAG 


AC P9 


GTCTACCAGGCATTCGCTTCATXXXXXGATGCCGCTG 


ACP 10 


GTCTACCAGGCATTCGCTTCATXXXXXTGGTCGTGCC 


ACP11 


GTCTACCAGGCATTCGCTTCATXXXXXCTGCAGGACC 


ACP 12 


GTCTACCAGGCATTCGCTTCATXXXXXACCGTGGACG 


ACP13 


GTCTACCAGGCATTCGCTTCATXXXXXGCTTCACCGC 


ACP14 


GTCTACCAGGCATTCGCTTCATXXXXXGCAAGTCGGC 


ACP15 


GTCTACCAGGCATTCGCTTCATXXXXXCCACCGTGTG 


ACP 16 


GTCTACCAGGCATTCGCTTCATXXXXXGTCGACGGTG 


ACP 17 


GTCTACCAGGCATTCGCTTCATXXXXXCAAGCCCACG 


ACP1 18 


GTCTACCAGGCATTCGCTTCATXXXXXCGGAGCATCC 


ACP 19 


GTCTACCAGGCATTCGCTTCATXXXXXCTCTGCGAGC 


ACP2D 


GTCTACCAGGCATTCGCTTCATXXXXXGACGTTGGCG 



found in the mirror carp (Figure 2). The full-length CcSLC 
cDNA consists of 1697 nucleotides with an 897 base pair 
open reading frame (ORF). The ORF encodes a peptide of 
298 amino acids with a molecular mass of 31.5 kDa and a 
theoretical isoelectric point of 7.49. There are seven amino 
acid positions that are different between the isoforms 
(Figure 2). SLC25a5 has six putative transmembrane do- 
mains, which are composed of three repeats of mitochon- 
drial carrier protein features by searching the Pfam 
database (Finn et al. 2010). SLC25a5 also contains one 
substrate binding site (Figure 2) and encodes an adenine 
nucleotide translocator (ANT2). ANT2 (also named as 
SLC25a5) is the most abundant mitochondrial protein 
(Itoi et al. 2005). 

Using the CcSLC2Sa5 mRNA sequence as a query, we 
searched our assembled genome database of Cyprinus 
carpio (not published), two copies of CcSLC2Sa5 exist on 
one scaffold (scaffold ID: 000000476), they are adjacent 
to each other and have the converse gene directions. 
By comparing the genomic and transcript sequences, the 
CcSLC25a5 has four exons and three introns (Figure 3). 
We then compared the carp gene to known zebrafish and 
human isoforms. Zebrafish have ANT1 {SLC25a4) and 
ANT2 (SLC25aS) genes, while humans have three ANT 
genes: ANT1, ANT2, and ANT3. Similar to the carp, the 
zebrafish ANT1 and the human ANT1, ANT2, and 
ANT3 genes contain four exons (Figure 3). The excep- 
tion is the zebrafish ANT2 gene which contains three 
exons (Figure 3). 

Phylogenetic analysis of cyprinid fishes and mammalian 
SLC25a5 isoforms 

To analyze the phylogenetic relationship of ANTs in 
cyprinid fishes, a phylogenetic tree was constructed 
using neighbor-joining from deduced amino acid se- 
quences (Figure 4). The six carp isoforms were clustered 
into one branch. They are more closely related to human 
and bovine ANT3, and more distantly to mammalian 
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KL-2 



JL-2 DL5000 KL-2 JL-2 





Figure 1 Differential gene expression detected in the Genefishing assay using primers ACP28 and ACP29. The PCR products from the 
genefishing assay were separated by agarose gel electrophoresis. Gene products that are differentially expressed in the mirror carp (KL-2) and 
common carp (JL-2) are shown by green arrows. DL5000 indicate molecular markers. 
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Figure 2 Comparing of the amino acid sequences deduced from six carp ANT2 transcripts. Dotted lines under the sequences indicate the 
positions of six putative transmembrane helices (Klingenberg 1 989) labeled segments l-VI. Three repeats of the mitochondrial carrier protein 
signature (Bof et al. 1999) are shown by dotted boxes. A substrate binding site is highlighted with shaded boxes. Amino acids that differ between 
the transcripts are indicated in blue. 
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Zebrafish ANT2 
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Figure 3 The gene structure of the carp, zebrafish, and human ANT genes. Coding regions in the exons are represented by black boxes. 
Introns are represented by solid lines. 



ANT Is (rat, mouse, bovine, human) and ANT2s (mouse, 
bovine, rat, human). The deduced amino acid sequence 
of the CcSLC25a5 protein shares a high degree of hom- 
ology with other known vertebrate SLCs (82% identity 
with mouse, bovine, rat, human). 



Tissue expression of the CcSLC25a5 

To define the timing and patterns of CcSLC2Sa5 expres- 
sion during early embryonic development, we performed 
in situ hybridization throughout the earlier developmen- 
tal stages. The SLC2Sa5 gene was expressed weakly in 
embryos from 4.5 (1/3 epiboly stage, Figure 5A) to 6 
(shield stage, Figure 5B) hpf (hours post fertilization). 
However, a strong expression signal appeared at the 70% 
epiboly stage (10 hpf, Figure 5C). The gene expression 
signal was concentrated in the embryo and undetectable 
in the yolk. SLC25aS was continuously expressed begin- 
ning 10 hpf (6-somites stage) to 16 hpf (14-somites stage) 
(Figure 5D-F) and continuing until 48 hpf (long-pec stage, 
Figure 5G) when the signal gradually weakened. Interestingly, 
CcSLC25a5 was specifically expressed in the scale placods 
as the scale pattern was formed (Figure 5H). 



Discussion 

Scales of the teleost fish are important skin integument- 
ary appendages distributed over the body surface in de- 
fined patterns. Natural variation in scale patterns exists 
among the common carp species. Mirror carp and Jianli 
are two carp varieties that have distinct scale patterns as 
a result of breeding selection. Scale initiation and mor- 
phogenesis are very complicated biological processes. 
The functional analysis of FGFR1 by reverse genetics 
showed that the mutation of this gene can lead to re- 
duced scales and abnormal fins (Rohner et al. 2009). The 
scale is one of the important agricultural traits for fishes 
and play important roles in physiology, defense, and 
adaptation to new environments (Sire et al. 1997). Some 
aquatic biologists are interested in scale development 
and launched some works using the reverse genetics. 
However, a complete knowledge of scale development is 
limited due to the longer sexual maturation period and 
the larger genome size for many species in fishes. There- 
fore, the molecular mechanisms underlying the scale ini- 
tiation and pattern formation remain unknown. Using 
Genefishing, we were able to identify seven genes that 
were differentially expressed in Jianli and mirror carp 
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Figure 4 A phylogenetic tree of carp and mammalian ANT isoforms constructed using the neighbor-joining methodology. The 

evolutionary distance between amino acid substitutions was estimated using Kimura's two parameter-method (Kimura 1980). The scale along the 
bottom axis shows the evolutionary distance of amino acid substitutions per site. 
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Figure 5 The CcSLC25a5 expression pattern by whole-mount in situ. Panels A-G show the expression pattern of CcSLC25a5 genes in 
different developmental stages. (A) 4.5 haf (1/3 epiboly), (B) 6 haf (shield), (C) 10 haf, (D) 12 haf, (E) 16 haf, (F) 24 haf, and (G) 48 haf after 
fertilization, respectively. (H) SLC25a5 is expressed specifically in the regions where the scales were developing in the skin. (I) There is no 
hybridization signal detected using the sense probe as the negative control. 



skin tissues (Table 2). These genes may contribute to the 
morphogenesis of the integumentary appendages, espe- 
cially the scales or fins, the functional analysis of these 
seven genes will be performed in future studies. 

ANTs are the most abundant mitochondrial proteins 
and mediate the exchange of ADP and ATP across the 
mitochondrial membrane. They link ATP production in 
mitochondria to its functional utilization for energy re- 
quirements outside mitochondria (Klingenberg 1981 and 
1989). Humans are known to have three ANT genes. 
ANT1 is predominantly expressed in the heart and skel- 
etal muscle and ANT2 is expressed in kidney and liver. 
ANT3 is expressed ubiquitously but is highest in the kid- 
ney (Stepien et al. 1992). The tissue distributions of dif- 
ferent ANT isoforms seem to reflect their functional 
differences and associate with tissue-specific energy me- 
tabolism (Stepien et al. 1992). 

Given their tissue specific roles in humans, it was of 
interest to examine expression profiles of ANT isoforms 
in various tissues in the carp. To date, there are no re- 
ports about ANT expression in fish scales. Herein, we 
report that ANT2 is specifically expressed in the scale 
placods. This result indicates that ANT2 may have an 
important role in scale development and morphogenesis. 
It has been suggested that fish ANTs have an important 
role in energy production, which is associated with 



temperature adaptation in fish. Guderley and Johnston 
(1996) reported that the uptake of ADP in isolated mito- 
chondria from cold-acclimated sculpin was higher than 
that of their warm-acclimated counterparts. CcSLC25a5 
possibly involved in the energy production in scale de- 
velopment for fishes, so it should be expressed in the 
process of scale development. Additional studies will 
provide new insights into the regulation of energy me- 
tabolism during scale development in fish. 

We next examined the expression profiles of CcSLC25a5 
in various tissues of the Jianli during development. 
CcSLC25a5 was weakly expressed in early stages of em- 
bryo development (4.5 to 6 hpf, Figure 5A-B). However, 
strong expression signals were detected in embryos from 
10 to 24 hpf (Figure 5C-F). These results imply that 
CcSLC25a5 has roles in embryo and scale development. 
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